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Abstract 

We present a detailed field theoretic description of those collective degrees of 
freedom (CDF) which are relevant to study macroscopic quantum dynamics of a 
quasi-one-dimensional ferromagnetic domain wall. We apply spin coherent state 
path integral (SCSPI) in the proper discrete time formalism (a) to extract the rele- 
vant CDF's, namely, the center position and the chirality of the domain wall, which 
originate from the translation and the rotation invariances of the system in ques- 
tion, and (b) to derive effective action for the CDF's by elimination of environmental 
zero-modes with the help of the Faddeev- Popov technique. The resulting effective 
action turns out to be such that both the center position and the chirality can be 
formally described by boson coherent state path integral. However, this is only 
formal; there is a subtle departure from the latter. 

1 Introduction 

Recent so-called nanostrucure technology enables us to study low dimensional magnetism 
in mesoscopic magnets from the quantum mechanical point of view [1-3]. Among others, a 
magnetic domain wall has attracted much attention both theoretically and experimentally, 
because it is expected to exhibit macroscopic quantum phenomena [4-11]. As a theoretical 
technique to evaluate the quantum dynamics of the domain wall, the spin coherent state 
path integral in the continuous-time formalism [12] is frequently used. However, as noted 
by some workers [13,14], it has some fundamental difficulties, which have been recently 
discussed in detail [15]. Furthermore, it is liable to lead to confusion concerning the 
interpretation of the collective degrees of freedom as has been pointed out in Ref.[16]. 
Hence, as yet there is no microscopic theory of the quantum dynamics of the domain 
wall. In order to pave the way for such a theory, this paper presents a field theoretic 
description of collective degrees of freedom by use of the spin coherent state path integral 
in the proper discrete-time formalism. 
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2 Model 



We consider a ferromagnet consisting of a spin S of magnitude S at each site in a quasi- 
one-dimensional cubic crystal (a linear chain) of lattice constant a. The magnet is assumed 
to have an easy axis in the z directions. Accordingly, we adopt the Hamiltonian 

n l N L 

<i,j> j 

where the index % or j represents a lattice point, < i, j > denotes a nearest-neighbor pair, 
N L is the total number of lattice points, and J is the exchange coupling constant, and K 
is the longitudinal anisotropy constants; J, K are all positive. 

Since we are interested in those transition amplitudes which are appropriate to describe 
quantum mechanical motion of a domain wall, we introduce a spin-coherent state [17] at 
each site, which is suited for a vector picture of spin. We denote a state of the system as 

IO = l6,6,---,U>:=(g)|^), (2) 
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where is a spin coherent state at the site j. The transition amplitude between the 
initial state and the final state |£p) can be expressed as a spin- coherent- state path 
integral in the real discrete-time formalism by the standard procedure of the repeated use 
of the resolution of unity (see, e.g., Ref.15 on which the present notation is based): 

N-l n l / ' \ 

(t F \e- iAT/n \b) = jim / n^(£>U;W)exp ( > (3) 

n=l jf ^ ' 

where N = T/e, e is an infinitesimal time interval, n represents discrete time, and the 
integration measure is 

MU):= (1 + \Z\*) 2 2m ' l^T = ^ ' (4) 

The action <S[£*, £] consists of two parts, <S C [£*, £] and <S d [£*, £], which are to be called the 
canonical term and the dynamical term, respectively. We shall be interested in those spin 
configurations whose scale of spatial variation is much larger than the lattice constant a. 
Accordingly, we take the spatial continuum limit in the action: 

S[C,i\ := S c [C^]+S d [C^}, (5a) 
: &\rf] - gf / L/2 ^ln (l + r(*,n)g(*,n-l)) 2 ^ 
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TS d [C,i] = -7-J2 e —H(e(x,n),Z(x,n-l)), (5c) 
"> J-LI* a 

H(C(x),v(x)) := 



-L/2 

2j^r(^)^(x) 



-ll^-^(i-r(^)) 2 + i 



(5d) 



where L is the length of the linear chain and J = J a 2 . 



3 Method of Collective Degrees of Freedom 
3.1 Kink configuration 

We begin by finding a domain wall configuration. It is determined by one of the static 
solutions {£ s (x), ^ s (x)} of the action £[£*,£]. They satisfy the following equations: 
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where A 2 = JS/K(S — 1/2). An obvious solution is the "vacuum" solution representing 
the uniform configuration in which the spins are either all parallel or all anti-parallel to 
the z direction. The other solution is the "kink" solution representing a domain- wall 
configuration in which the spins at x ~ +oo are parallel to the z direction, the spins 
at x ~ — oo are anti-parallel to the z direction, and there is a transition region (i.e., a 
domain wall) of width A; 

f (x) = exp (^-—J-^ + i( l>oj > f 0) = ex P ^ — - «0o^ , (7) 

where Q and 0o are arbitrary real constants. Q is the center position of the domain wall, 
and 0o is a quantitative measure of the chirality of the domain wall with respect to the x 
axis (Fig. 1); the wall is maximally right-handed if 0o — 7r /2 and maximally left-handed if 
<po = — vr/2, while it has no chirality if 0o — 0. The range of 0o is chosen as — tt < <po < ti, 
with <pQ = tt and 0o = — tt representing the same situation. {£ s (— x), <£ s (— x)} is also a 
solution representing a domain-wall configuration. However, this as well as the vacuum 
solution belongs to a sector different from that of (|^). Since a transition between different 
sectors is forbidden [18], it is sufficient to consider only the sector (|7|) for the purpose of 
studying the dynamics of a domain wall. 



3.2 Collective degrees of freedom and environment 

Study of the domain-wall dynamics is facilitated by introducing relevant collective degrees 
of freedom. We note two kinds of invariance possessed by Eq. ([)]). One is the translation 
invariance in the x direction. The other is the rotation invariance around the z axis. These 
invariances are embodied by the arbitrariness in the choice of Q and 0o, respectively, in 
(0). Hence we elevate them to dynamical variables Q(n) and <fio{n) [8,10,18]. To deal with 
these two dynamical variables (collective degrees of freedom), it is convenient to define 



z(n) := q(n) + i(f>o(n), z*(n) := q(n) - i(f>o(n), 
q(n)=Q(n)/X, n = l,2,...,JV-l. 



(8) 




Figure 1: Domain walls with three chiralities (quoted from Ref.10); (a) right-handed wall 
(0o = 7r/2), (b) left-handed wall (0o = ~ 7r /2), and (c) wall with no chirality (0o — 0). 
Circles in (a) and (b) drawn to guide the eye lie in the yz plane, while the spins lie in 
the zx plane in (c). The quasi-one-dimensional direction of the crystal is here aligned 
with the spin hard axis for ease of visualization. A different alignment, which may be 
the case for a real magnet, does not affect the content of the text; for instance, one could 
rotate all the spins by it/ 2 around the y axis if the dominant anisotropy originates from 
the demagnetizing field. 

By use of these variables, original integration variables £(x,n) and £*(x, n) may be de- 
composed into the domain-wall configuration and the deviation from it: 

£(x,n) = i s {x;z{n))+fj{x,n]{z}l), (9a) 

C(x,n) = e(x;z*(n))+fj*(x,n;{z*}l), (9b) 

where 

C{x;z{n)) :=exp{-x/\ + z{n)) , f (x; z*{n)) := exp {-x/X + z*(n)) , (10) 

and we use the notation;{2:}^ = {-2*}™ := (z*(n), z(m)). At both ends of the discrete 
time (n = or n = N), we define 

z(0) = zi := qi + i(p h z(N) = z F := q F + i(f) F , (11) 
V (x,0;{z}° o )=r 1 *(x,N;{z}%) = 0, (12) 

where z\ represents the center position q\ and the chirality (pi of the domain wall in the 
initial state, and z F those in the final state. 

The variable fj, which are to be called environment around the domain wall, could be 
expanded by a set of some mode function. The expansion is expected to contain the 
zero-modes, which originate from the translation and the rotation invariances and should 
be eventually eliminated in order to avoid overcounting the degrees of freedom. For this 
reason, we expand the environment with respect to a set of mode functions {ipk} as 



ri(x,n;{z}™) = r] (ri)ip (x; z(n)) + r)(x, n; {z}%), (13a) 

fj*(x,n;{z*}Z) = V * (n)r (x;z*(n))+r]*(x,n;{z*}l), (13b) 

V (x,n;{z} n n ) = J2'vk(n)il; k fr{z}n), (13c) 
k 

V *(x,n;{zTn) = E'^( n )^M**}n)> ( 13d ) 
k 



where Yl'k denotes summation over the modes excluding the zero-modes. The zero-modes 
function is proportional to the first derivative of £ s (x; z(n)) with respect to x; 

i , / w ^ d^{x\ z(n)) /-, ,x 

z{n)) = —AX y ^ JJ , (14) 

where A is a real normalization constant. To confirm this relation, one may consider the 
sum of the domain wall configuration and the zero- mode part: 

C(x; z(n)) + i] Q (n)tjj (x; z{n)) 

/ n, / n d£ s (x: z(n)) 

= C{x; z(n)) - AXr) (x) v ^ y 11 

— C( x — A\rj (n); z(n)) 

= exp \~ + q{n) + AUr] {n) + i(M n ) + 49tyo(n))l . (15) 
L A 

Hence, it is clear that Eq.([L4]) gives the zero-mode eigen function, and that the real and 
the imaginary parts of r] (n) are the zero- modes corresponding to the translation and the 
rotation modes, respectively. We choose {ipk} so that {tp Q , ip^} forms a orthonormal set; 

L/2 dx 

—f{x- {z}l)r (x; z\n))%{x- z{n)) = 1, (16a) 

-L/2 a 

L/2 dx 

—f{x; {z} n n Wo(x; z*(n))Mx; {z}D = 0, (16b) 

-L/2 a 

L/2 dx 

— /fa U}rM(x; {z*Yn)Mx\ UVn) = <W, (16C) 
-L/2 a 



f(x; {z}n){M x \ z{n))%{x\ z*(n)) 

+£V*te {zTj] = S {^-) ■ (16d) 

where f(x; {z}™) is a real weight function to be fixed later. This weight function neglects 
a non-linear character of the coherent state path integral. 

3.3 Faddeev-Popov type identity 

In order to introduce the collective degrees of freedom and eliminate the zero modes, we 
adopt the Faddeev-Popov method to the spin-coherent state path integral. To do so, we 
invoke the following Faddeev-Popov type identity: 

dzW^ (1?a) 



where 



dz{n)dz*{n) 

— = dq(n)d(po{n), (17b) 



and AK'£*l({z}£) is the Faddeev-Popov type determinant 

d($ ] ({z} n n ),g* Q r \{zTn)) 



d(z(n), z*(n)) 



(1? 



with Rg^'^ \{z}™) and Ig^^ J ({^}^) being the real and the imaginary parts of fl^' ({#}?■ 
defined as 



-/(x; {*}£){£(*, n) - f (x; «(n))}^(x; **(n)), (19a) 

-L/2 a 

^({^JS := tf (MS)}*- (19b) 



3.4 Transition amplitude 

Hereafter we consider transition amplitudes between domain-wall states. Namely we take 



\^) = h):=6Q\e(ja;z p )), /3 = I,F. 



(20) 



in Eq. @. Inserting the identity (|17a|) into the right-hand side of (H) at each discrete 
time, we get 
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AT-l 
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xex P ( -,s[r,e] 



(21) 



At this stage, we substitute (§) with (|16|) for £ and £*. Then (fl9| ) may be replaced by 

g [ f +ii] ({z} n n) = Vo(n), gf +r] ({z*} n n) = %»• ( 22 ) 

Thus, 



^0 



({«}») = ^o(n), j a 



({*}£) = ^o(n), 



(23) 



which are zero-mode degrees of freedom. The integration measure can be rewritten as 

d^(x,n)dC{x,n) 2S + 1 
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fc=0 



2m 

di] k (n)dr]l(n) 
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(24) 



where 



[x,n;z,fj] :-- 



{1 + (e(x; z*(n)) + fj*(x; {z*}»))(?>(x; z(n)) + fj(x; {z}%))} 2 
Thus, the transition amplitude (|2T| ) reduced to 
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(26) 



With the help of the delta function, we can immediately integrate out the zero-modes at 
each discrete time. Furthermore, in the case of large spin S ^> 1, The factors contributing 
to the integration measure takes the simple form [19]: 



Z[x,n;z,r)} = 1 + 0(S- 1/2 ) 



(! + (S- 1 ' 2 )) 
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In this way we find 

(^le-^lzr) 



JV-l 



lim / I f 

N^oo I - LJ - 



n=l 



dr) k {n)dr)* k {n) dz(n)dz*(n) 



2m 



2m 



xex P (-5[f + v*,C + V\ 



:28) 



This formula determines the quantum dynamics of the domain wall in the path integral. 



3.5 Effective Action 

The next step is to expand the action 5[£ s + 77*, £ s + 77] with respect to the environment. 
The expanded action consists of three parts. The first part <S c [z*,z] contains only the 
collective degrees of freedom. This part was estimated in detail in Ref. [16]. The second 
part S c [f]*,r]] contains only the environment, and the third part S c - e [z*, z, 77*, 77] represents 
the interaction between the collective degrees of freedom and the environment; 

S[? + V*, e + v] ■= S c [z*, z] + S c [ V *,ri} + S c . e [z*, z, ?f, rj\. (29) 

Now, we proceed to find on equation to determine the environmental mode functions 
{ipk}- Such an equation may be obtained from the dynamical term whose collective part 
gave the static equation(|6|). Accordingly, we examine the dynamical term of «S e [77*, 77] , 



which we denote by [r]*,rj\: 
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where 



~n _ ~*m 
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x z*{n) + z(m) 
A 2 



9{x\{z} n n ) 



7y(x,n (30a) 



(30b) 



1 + e~ 2 ^ 



Terms with proportional to 77*77* or rji] vanish due to the static equation^ . Note that 
this cancellation is not a generally guaranteed theorem. It turns out to be conevenient to 
write {ipk,4>k} m ( P-3cp and ( |13d[ ) as follows: 



i/) k (x;{z}$) := 2e- x/x+z{n) cosh 
Putting ( TO , ( |T3c| ) and (gT^) into (gTFj), we find 
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(31b) 



V / n =l fc'.fc 



X / — Vfc'(^) 
-L/2 a 



<9 2 



1 



dx 2 X 2 



+ — M-2sech" - 



x 



Pk{X) 



(32) 



where we replaces x/X — (z*(n) + z{n))/2 by x/X because L/X ^> 1. Thus, we are led to 
demand <pk(x) to obey the eigenvalue equation 



(33) 



where u;^ is the eigenvalue. This solution to this equation is well known in the literature; 

(p k (x) = N k {-ikX + tanh(x/A)) e ikx } (34a) 

where 



(k 2 X 2 + 1) 



-1/2 



tu k = k 2 + A~ 



(34b) 



These eigenfunctions giagonalize S^[r}*,rj\ as 



N 



(35a) 



n=l k 



where 



(35b) 



The last equation is the dispersion relation of the environmental eigenmode. Note that 
ujk is not zero in the limit of k — » 0. 

What remains is to decide the weight function f(x; {z}™). With (|3T|), the condition 
(16) takes the form 

rL/2 j 

-f(x, {zT n )Ae- 2 ~< cosh 2 {xl) vl(x; {^1^ = (36) 

-L/2 a 

This is satisfied if we choose 

1 



(1 _|_ e -2x/\+z*(n)+z(n)y 



(37) 



It is easy to check that (|16b|) is also satisfied by this choice. Accordingly, the normalization 
constant A of the zero- mode function ipo{x] z (n)) as given by flT4] ) may be determined by 
use of (|37|) and ( |16a| ): 



l = A 2 



L/2 ^ x e -2x/\+z*(n)+z(n) 
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Thus, 
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(38) 



(39) 



The closure ( |16d|) can also be shown to hold with these choices. 

The remainder of the action can be calculated in the same manner [19], and we obtain 
the effective action in the following form: 
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(40c) 
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where Az(n) = z(n) — z(n — 1), and 



A f L/2X sinhf , x , x 

2 « J-L/2A COSh (z + § J 

which is a function of the collective degrees of freedom, represents a non-linear interaction 
with the environment. It is seen that the above form of the effective action is formally the 
same as that obtained in a boson coherent state path integral. Hence, it is concluded that 
the quantum dynamics of the domain wall is formally represented by that of a "boson" 
z interacting with environmental bosons {77}. However, it is to be remenbered that the 
imaginary part of the "boson" z is an angular variable. This circumstance can cause a 
subtle departure from the case of a boson. Details of these subtleties as well as a concrete 
evaluation of transition amplitudes are left for a future work. 
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